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We investigate the capabilities of a quantum computer based on cold trapped ions in presence 
of non-dissipative decoherence. The latter is accounted by using the evolution time as a random 
variable and then averaging on a properly defined probability distribution. Severe bounds on com- 
putational performances are found. 
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Q I. INTRODUCTION 

<N \ 

A quantum computer processes quantum information that is stored in quantum bits (qubits) Q. If a small set 
of fundamental operations, or "universal quantum logic gates" can be performed on the qubits, then a quantum 
. computer can be programmed to solve an arbitrary problem || . Essentially, the quantum computation can be viewed 
(f) ' as a coherent superposition of digital computation proceeding in parallel. The explosion of interest in quantum 
t-H , computation can be traced to Shor's demonstration that a quantum computer could efficiently factorizc large integers 



Cirac and Zoller Q have made one of the most promising proposal for the implementation of a quantum computer. 
A number of identical atoms are stored and laser cooled in a linear radio-frequency quadrupole trap to form a quantum 
, register. The radio-frequency trap potential gives strong confinement of the ions in the directions transverse to the trap 
axis, while an electrostatic potential forces the ions to oscillate in an effective harmonic potential in the axial direction. 
After laser cooling the ions become localized along the trap axis with a spacing determined by their Coulomb repulsion 
and the confining axial potential. The normal mode of the ions' collective oscillations which has the lowest frequency, 
is the axial center of mass (CM) mode, in which all the trapped ions oscillate together. A qubit is the electronic 
ground state \g) and a long-lived excited state |e) of the trapped ions. The electronic configuration of individual ions, 
and the quantum state of their collective CM vibrations can be manipulated by coherent interactions of the ion with 
a laser beam, in a standing wave configuration, which can be pointed at any of the ions. The CM mode of the axial 
vibrations may then be used as a "bus" to implement the quantum logical gates. Once the quantum computation has 
been completed, the readout is performed through the mechanism of quantum jumps. Several features of this scheme 
^ , have been demonstrated experimentally, mostly using a single trapped ion Q . 

The implementation of a large-scale quantum computer is recognized to be a technological challenge of unprecedent 
, proportions. In fact the qubits must be easily manipulated, but must be also well isolated from decohering influence 
which places considerable limitations on the capabilities of such devices ||. Practically, there are two fundamen- 
tally different types of decoherence during a computation: the dissipative one, due to the spontaneous decay of the 
metastable states |e), and the non-dissipative one due to random phase fluctuations of various nature. While the first 
has been considered carefully , by also developing strategies to reduce spontaneous decay, like watchdog stabilization 
H , the second does not received much attention Ref. |j| . 

On the other hand, non-dissipative decoherence seems actually dominant in trapped ion based experiments p0[ , 
whose results have been recently well explained , both qualitatively and quantitatively, by simply considering the 
time as a statistical variable and then averaging over a properly denned probability distribution ]l2| . The developed 
approach results model-independent, hence quite ductile. Then, our aim here, is to study the performances of an 
ion-trap quantum information processor by using this approach. Practically, the latter allows us to account for the 
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fluctuations of the Raman laser intensity |13|. Instead, in [|9| a specific model is used, i.e. a random phase drift is 
added by hand at the end of each qubit rotation while in our case the stochasticity is already included in the time 
evolution. Moreover, our analysis is based on exact parameters values fitting experimental data pd[ . 

The paper is organized as follows. In Section II we shall review a general model for non-dissipative decoherence. 
In Section III we shall apply this theory to fundamental logic operations. In Section IV we evaluate the limits of a 
quantum computer due to this type of decoherence. Finally, Section V is for concluding remarks. 
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II. A GENERAL FORMALISM FOR NON-DISSIPATIVE DECOHERENCE 



In this section we review the formalism describing non-dissipative decoherence derived in Refs. |12|]. It is based 
on the idea that time is a random variable or, alternatively, that the system Hamiltonian (therefore its eigenvalues) 
fluctuates. This leads to random phases in the energy eigenstates representation. Then, the resulting evolution of 
the system must be averaged on a suitable probability distribution, and this leads to the decay of the off-diagonal 
elements of the density operator. 

Let us consider an initial state p(0) and consider the case of a random evolution time. Then, the evolved state will 
be averaged over a probability distribution P(t,t'), i.e. 

p(t)= / dt' P(t,t') P (t') , (1) 
Jo 

where p(t') — exp{— iLt'}p(0) is the usual solution of the Liouville-Von Neumann equation with L . . . = [H, . . .]/%, 
One can write as well 



p(t) = V L (t)p(0), (2) 

where 



V L (t) = / dt'P(t,t')e- tLt . (3) 
Jo 

In Ref. ]l2| ], the function P(t,t') has been determined to satisfy the following conditions: i) p(t) must be a density 
operator, i.e. it must be self-adjoint, positive-definite, and with unit-trace. This leads to the condition that P(t,t') 
must be non-negative and normalized, i.e. a probability density in t', so that Eq.(Q) is a completely positive map; ii) 
Vj,(t) satisfies the semigroup property Vl(£i + 1 2 ) — Vl(£i)Vl(£2), with t\,t% > 0. These requirements are satisfied by 



and 
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where the parameter r naturally appears as a scaling time. The expression (^|) is the so-called T-distribution function, 
well known in line theory. Its interpretation is particularly simple when i/r = k with k integer; in that case P(k,t') 
gives the probability density that the waiting time for k independent events is t' and r is the mean time interval 
between two events. Generally, the meaning of the parameter r can be understood by considering the mean of the 
evolution time (t') = t, and its variance (t' 2 } — (t 1 ) 2 = rt. 

When t — > 0, P(t,t') — > S(t — t') so that ~p(t) = p(t) and Vl(£) = exp{— iLi} is the usual unitary evolution. 
However, for finite r, the evolution operator Vi(t) describes a decay of the off diagonal matrix elements in the energy 
representation, whereas the diagonal matrix elements remain constants, i.e. the energy is still a constant of motion. 
In fact, in the energy eigenbasis, Eqs.(||) and ([!]) yield 

p n , m (t) = p n>m (0) , (6) 

where 

7«,m = log (l + c4, m r) , (7) 

v n m = - arctan (ui n m r) , (8) 

T 

with Tiuj n m — {E n — E m ) the energy difference. One can recognize in Eq.(^), beside the exponential decay, a frequency 
shift of every oscillating term. 

The phase diffusion aspect of the present approach can also be seen in the evolution equation for the averaged 
density matrix p(t). Indeed, by differentiating with respect to time Eq.(|^) and using (||) one gets the following master 
equation for ~p(t) 
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j t p{t) = -~\og{l + iLT)T>{t)- (9) 
It is worth noting that by expanding the logarithm up to second order in t, one obtains 

j-p{t) = - l - [H,p(t)} ~ ^2 [B, [H,p(t)}} , (10) 

which is the well known phase-destroying master equation ||l4| . Hence, Eq. (||) appears as a generalized phase destroying 
master equation taking into account higher order terms in r. Nonetheless, the present approach is different from the 
usual master equation approach, in the sense that it is model independent, non perturbative and without specific 
statistical assumptions. 

The theory well describes also Hamiltonian fluctuations. Let us consider, for instance, a Hamiltonian H (t) — fl(t)H, 
where tt(t) is a fluctuating parameter with mean value O. Whenever H = H/Q induces oscillations, we can accounts 
for the fluctuations of the frequency O(i) by writing |ll| 

p{i) =exp{-iL4(t)}p(0), =S> p(t)= J dAP(t,A)e- llA p(0), (11) 
where L = [H,.. .]/%, and A(t) = /J d£n(£) is a positive dimcnsionless random variable. Furthermore, 

P(M)4^T^ ■ (12) 



fir r(t/r) \CIt / 

The first two moments of P(t, A) determine the properties of the fluctuating frequency £l(t), 

(A) = fii , a 2 (A) = (A 2 ) - (A) 2 = fl 2 tT , (13) 

that is, the frequency f2(f) is a white, non-gaussian (due to the non-gaussian form of P(t,A)) stochastic process. In 
fact, the semigroup assumption made, implies a Markovian treatment in which the spectrum of the fluctuations is 
flat in the relevant frequency range. This in particular implies that we are neglecting the dynamics at small times, 
of the order of the correlation time of the fluctuations. The non-gaussian character of P(t, A) can be traced back to 
the fact that P(t,A) must be definite and normalized in the interval < A < oo and not in — oo < A < oo. This is 
another fundamental difference with respect to Ref . |l , where a gaussian character for the random phase drift errors 
is assumed, while we make no a priori statistical assumptions. Indeed, the properties of the probability distribution 
P are derived only from the semigroup condition, and it is interesting to note that this condition yields a gaussian 
probability distribution as a limiting case. In fact, from Eq.(||) one can see that P(t,t') tends to become a gaussian 
in the large time limit f>r. 

In the case of Rabi oscillations, A becomes proportional to the pulse area, while r gives an estimate of the pulse 
area fluctuations, since it corresponds to a fractional error of the pulse area J a 2 (A) j (A) = \frjt. 



III. CHARACTERIZATION OF A QUANTUM PROCESS 



Quantum computation ideally corresponds to a physical process I'J'out) = U\$\r) where a given input state is 
mapped to an output state by a unitary transformation U. It has been shown that any quantum computation can be 
decomposed into one-bit gates and a universal two-bit gate which involves entanglement operation on two qubits [fl5[ . 

In order to characterize a physical process in a quantum system, we follow Ref. flo] . The aim is to characterize the 
process given as a "black box" , by a sequence of measurements in such a way that it is possible to predict what the 
output state will be for any input state. 

We assume that the system is initially prepared in a pure state 

N 

|*5n)=X) C *l i )' Pin = l*in)(*in|, (14) 

i=0 

where |0), |1) . . . , \N) are orthogonal states spanning the Hilbcrt space of allowed input states with dimension N + 1. 
Then, in the ideal case, the output state can be written as 
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|*out)(*out| = Pout = ^2 fyC*, Ri< :i , (15) 
i,i'=0 

where 

Ri',i = U\i)(i'\U\ (16) 

are system operators not depending on the input state. 

Whenever the unitary transformation is given by the time evolution operator U(t), the initial state p(0) represents 
the input state while the evolved state p(t) represents the output one. Since in Section II p(t) is replaced by ~p{t), 
here we have to replace Ri\i(t) by 

poo 

Ri>,i(t)= dt'P(t,t')U(t')\i){i'\U\t'), (17) 



Now, in order to see to what extent the real physical process approaches the ideal one, we use as parameter the 
fidelity 

F=Tr{p(t)p(t)} am , (18) 

where the subscript "ave" indicates the average overall possible input states. Obviously a gate fidelity close to one 
indicates that the gate was carried out almost ideally. 

In what follows we will see that the fidelity T ([18]) can be expressed in terms of the tensor 



Fp^ifltfR^Ulj). (19) 



A. One-bit gate 

In the case of one bit N = 1, and the vector basis is given by {\i = 0) = \g), \i = 1) = |e)}. The individual rotation 
acting on a single ion can be performed using a laser frequency on resonance with the internal transition and with the 
equilibrium position of the ion coinciding with the antinode of the laser standing wave. The Hamiltonian describing 
such a interaction is 

H=^[\e)(g\e-^ + \g)(e\e^] , (20) 

where f2 is the Rabi frequency and 4> the laser phase. Eq.(p0|) leads to the following evolution 

U(t)\g) = cos [Vtt/2] \g) - ie** sin [Q.t/2] |e) , (21) 
U(t)\e) = cos [Vtt/2] |e) - ie'^ sin [flt/2] \g) , (22) 

which effectively corresponds to a single qubit rotation. It is clear from the above evolution equations that fluctuations 
in the Rabi frequency can be accounted by simply replacing with A(t) — J dt;tt(t;) (where now Q indicates the 
mean Rabi frequency) accordingly to the arguments of Section II. Then, similarly to Eq.(fi~7]), it is possible to evaluate 



Ri>,i(t)= / dAP(t,A)U(A)\i)(i'\U^(A), (23) 
Jo 

and the tensor elements ( p^ ) (see Appendix I) . In terms of these elements the fidelity results 

^lE^+i E ( F ii +F ii) ■ (24) 



i=0 i^j=0 

In this case the fidelity depends on time which practically determines the amount of rotation, however, in the next 
Section we shall consider a 7r rotation, i.e. t = n/Q. 
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B. Two-bit gate 



In the case of two-bit N = 3, and the vector basis is given by {\i = 0) = |<?)i|g)2 , \i = 1) = |.9)i|e)2, |i = 2) = 
|e)i|g)2, |* = 3) = |e) i |e)2 where labels 1, 2 indicate the two bits. Beside that also the vibrational ground state |0) 
is employed. 

In this case the laser frequency is chosen tuned to the first motional sideband and the equilibrium position of the 
ion coincides with the node of the laser standing wave. The Hamiltonian describing such interaction is Q 



f [|e) n ( ff |ae-^ + |<?)„(e|at e ^] , q = . 



n' = 



^[\e') n (g\ae-^ + \g) n (e'\ah^] 



7k 



(25) 



Here, and a are the creation and annihilation operators of the CM phonons, Q is the Rabi frequency, 6 is the laser 
phase, and rj is the Lamb-Dicke parameter. The index q = 0, 1 refers to the transitions excited by the laser, |<?) <-> |e) 
or \g) <-> |e'), which depend on the laser polarization. Instead the index n refers to the n-th ion on the trap (1, 2 in 
our case). Moreover, the factor \/N a , where N a indicates the number of trapped ions, appears as a consequence of 
the Mossbauer effect B. 

The universal two-bit gate, defined by |l5|] 



|ei)i|e 2 )2^(-l) eie2 |ei)i|e2)2, (ci,2 = 0, 1) 
can be realized in three steps by means of Eq. (p5| ) : 

• A 7r laser pulse with polarization q = and = excites e.g. the first ion. The evolution will be 



(26) 



exp 



(|e)i<<?|a+l.9)i(e|a f ) 



(27) 



with t\ = 7r/ST. 



• The laser directed on the second ion is then turned on for a time of a 2n pulse with polarization q = 1 and 
6 = 0. The evolution will be 



Un{t 2 ) = exp 



-i^h(\e') 2 {g\a+\gh{e'\J) 



(28) 



with t 2 = 2ir/n'. 

• A 7r laser pulse with polarization q = and = excites again the first ion. The evolution will be 



exp 



Vt' 

-*Y*3(|e)i<s|a + | 5 )i(e|aT) 



(29) 



with in 



tt/Q'. 



Also in this case to account for the fluctuations in the Rabi frequency, we can introduce a stochastic variable A' (t) = 
J Q d£ f2'(£) = J Q d£ I \J N a (where now Vt! is related to the mean Rabi frequency ft as in (^5|)). Then, we have 

/>oo poo poo 

R\ H = dA' 3 dA' 2 dA , 1 P{h,A'z)P{t 2 ,A' 2 )P{t l ,A' l ) 

Jo Jo Jo 



which are used to calculate the tensor elements (|l^) (see Appendix II). Finally, the fidelity results 

3 , 3 



IsTf" + — V (f*. 



p3 l 

l 3 



(30) 



(31) 



i=0 

In this case T does not depend on time since the gate operation is realized with the above definite steps. 
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IV. BOUNDS FOR QUANTUM COMPUTATION 



If one accounts the recent breakthroughs in the real of fault tolerant quantum computation [|l7| , an arbitrarily large 
quantum computations can be performed accurately provided to have an high degree of accuracy on the single gate 
operation. This sets an obvious figure of merit for quantum computation technology, namely, the expected probability 
of error in one quantum gate, which should be of the order of 10~ 6 JTs| . 

In Fig.l we show the deviation from the perfect gate fidelity as a function of the fractional error of the pulse area. 
In particular the lower straight line concerns the one-bit gate for a ir rotation (pl|), and the upper straight line the 
universal two-bit gate (|3l|). It is known that the value of 10 -8 s gives the best fit for the experimental data 
of Ref. |ic|| , i.e. £1 « 10 5 s _1 , r/ « 10 -1 . It means to have Qt m 10~ 3 whose corresponding value of accuracy, i.e. 
1 — J 7 , is quite far from the desired one, i.e. 10 -6 . In the case of universal two-bit gate, we have a better value of 
accuracy (upper straight line) already for Q't slightly less than fir. It turns out that the main limitations rely on 
the single bit rotation instead on universal two-bit gate. This fact, though counterintuitive, can be easily understood 
if one consider the expression ( |l3| ) for the two cases. It results a 2 (A) /a 2 (A') — \JN a /rj > 1 for 7r-pulses. This means 
that the coupling of the qubits with the vibrational degree of freedom makes the two-bit gate less affected by the 
noise with respect to the one-bit gate. On the other hand, the latter is built up with on resonance pulse, hence it is 
more sensitive to noisy effects, differently from the two-bit gate where off resonant pulses are used. However, since 
actually r]/-\/N a « 10 _1 , also in the two-bit gate the accuracy falls very far from the desirable value. 

Based on these results, we can state that quantum information processing on a large scale is unrealistic. As a 
matter of fact, the fault tolerant quantum computation requires a value fir w 10 -6 , i.e. r « 1CP 11 s for il w 10 5 s _1 . 
Since the finite value of r is related to the Rabi frequency fluctuations (or in turn to laser intensity fluctuations) [[l~l| , 
this means to improve the laser stability by a factor 10 2 at least! 

It is also to remark that, within the presented non-dissipative decoherence theory, the linear behavior of the gate 
fidelity (Fig.l) is typical of the limit fir -C 1. 

A rough estimation of the capabilities of a quantum computer in the presence of non-dissipative decoherence can be 
also made with the following arguments. Let us consider the optimistic case of a single run of the Shor's algorithm, 
and the use of the universal two-bit gate operations. Then, the non-dissipative decoherence theory |l2| shows that in 
the limit of Q't <C 1 the decay rate can be written as 7 = 2W 2 t (see e.g. Eq.(m)). Let us now suppose to factorize 
a L-bit number. Then, the number of trapped ions N a should be at least 5L M. The time needed for a single run 
of the Shor's algorithm is given by the time required for an elementary logical operation Air^/bL/rifl multiplied by 
the required number of elementary operations, approximately (10 L) 3 m. Of course this product should be much 
less then the docoherence time 7 -1 . Nevertheless, by using the actual experimental values for parameters (and the 
corresponding value of r « 10 -8 ), it results clear that the factorization of even a four-bit number results impossible. 

V. CONCLUSION 

In conclusion we have studied the limitations imposed by non-dissipative decoherence on quantum computation. 
Practically, we have seen that non-dissipative decoherence actually constitute a serious impediment to realize quantum 
computer beside dissipative decoherence. We have used a model able to accurately describe the decoherence phenom- 
ena on ion trap based experiments caused by the fluctuations of classical quantities. We have shown that large scale 
computation seems impossible with the present proposals. Our results indicate that even a computationally modest 
goal will be extremely challenging experimentally jljj . 

Although the conclusions of this paper are rather pessimistic with regard to the practical application of quantum 
computers for actual computation, there are applications requiring much fewer operations which are worth considering 
PI- 
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APPENDIX I 



In this Appendix we explicitly calculate the tensor elements ( |l9| ) for one-bit gate which are useful to calculate the 
gate fidelity. We make use of Eqs.(0), @ and (|l]), (§|) obtaining 



F ° ° = Fll = F™ = FU = / dA P(t, A) cos 2 [(A - fit)/2] 

Jo 



i arctan(£lr) 



+ sin(f2i) sin 



t arctan(£lr) 



2(1 + ^ 2 7 
= (1 

The obtained results are valid in any regime even if we are interested in f2r <C 1 and t — tt/Q. 



pii _ F oo 
r oo — r ll 



(1 - 1$S) 



(32) 
(33) 



APPENDIX II 



In this Appendix we explicitly calculate the tensor elements (|19|) for two-bit gate which arc useful to calculate the 
gate fidelity. Let us first examine the effects of evolution operators (|27j), (|28[), (|29| ) on the two-bit vector basis. For 
i = 0, we have 



For i = 1, we have 



U III U II U I \g) 1 \g) 2 \0) = \g) 1 \g) 2 \0) 
UmUuUi^ehP) = |fl)i|e) a |0). 



For i = 2, we have 

C////t///t/ / |e> 1 | 5 > 2 |0) =cos [fl'ti/2] {cos [fl't 3 /2] |e)i|ff>2|0) - i sin [n't 3 /2] |g)i| 5 ) 2 |l)} 

-zsin [n't x /2] cos [n't a /2] {cos [n'* 3 /2] |g)i| 5 ) 2 |l) - isin [n'* 3 /2] |e)i|g) 2 |0)} 
- sin [fl'ti/2] sin [0'i 2 /2] |fl)i|e')a|0) . 

For i = 3, we have 

t/ /J/ ^ J [/ / |e) 1 |e) 2 |0)=cos[fi / ii/2]{cos[n%/2]|e) 1 |e) 2 |0)-isrn[n'f 3 /2]| 5 ) 1 |e) 2 |l)} 

-i sin [fl% /2] {cos [n't 3 /2] | g) i |e> 2 1 1) - t sin [fi / t 3 /2] |e) x \e) 2 10) } . 
Then, by using the above results in Eqs.(|l7j) and ([ji]), we get 

TTlll _ pio _ pOl _ -1 



poo 

?2 2 



F|| = [C 2 (7r/Q')] + [S 2 (7r/fi')] C2(27r/n') " 2 [Z(7r/n')] Ci(27r/I2') 



K 3 3 _ ^ ,_/rVM 2 



33 - [c 2 (n/n')Y + [s 2 (n/n')} 2 2 [z(V^')] z , 

02° = *?Z = F 2 02 = F\l = [Cifr/fi')] 2 - [^(vr/rOf C 1 (27r/fi') . 



p^U _ F Zi _ pUZ _ p 

= Ffi = f° 3 = f 3 \ 3 = - [c^/n't + [^(Tr/n')]" , 

F 2 3 3 2 = FH = - [C 2 (ir/n')] 2 - [S 2 (n/n')] 2 d^/fi') + ^(V^Of + [W^Of C^Tr/O') , 

p0 _ p0 _ p0 _ ell _ 771II _ pill _ p2 2 _ p2 2 _ p2 2 p3 3 
^11 -^22 -^33 —-^00 -^22 ~^33 — — i*i l — ^ — 



p^ ^ p2 ^ p3 3 p3 3 p3 3 n 

■^11 —^33 —-^00 —^11 —-^23 — U ' 

t arctan(ri'r) 



where 



C 2 (t) 



[ dA'P{t, A') = 2(1 + J V2)t/2T {(1 + tfV)'/* + 



/>OC 

Z(t) = / dA' P(t, A') sin[A'/2] cos[A'/2] = 
Jo 

/>oc 

S 2 {t) = / dA' P(t, A') sin 2 [A'/ 2] = 1 - C 2 (i) , 
Jo 



1 



2(1 + 0' 2 T 2 )*/2t 



sm 



t arctan(f2V) 



d(t) 



dA' P(t, A') cos[A'/2] 
dA' A') sin[A72] 



1 



(l + fi' 2 T 2 /4)*/2r 
1 

(l+fi' 2 T 2 /4)*/2r 



t arctan(fi'r/2) 



sm 



t arctan(i7'r/2) 



(34) 
(35) 



(36) 

(37) 

(38) 
(39) 
(40) 
(41) 
(42) 
(43) 
(44) 

(45) 
(46) 
(47) 
(48) 
(49) 
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The obtained results are valid in any regime even if we are interested in O'r <C 1. 
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FIG. 1. The quantity 1 — T vs the fractional error of the pulse area is shown in a log-log plot. The lower straight line 
concerns a one- bit gate (after a 7r rotation). The upper straight line concerns the universal two-bit gate. 
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